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Abstract 

The quantum "Zeno" time of the 2P-1S transition of the hydrogen atom is 
computed and found to be approximately 3.59-10~ 15 s (the lifetime is approximately 
1.595 • 10 _9 s). The temporal behavior of this system is analyzed in a purely 
quantum field theoretical framework and is compared to the exponential decay 
law. 



Unstable systems decay according to an exponential law. Such a law has 
been experimentally verified with very high accuracy on many quantum me- 
chanical systems. Yet, its logical status is both subtle and delicate, because 
the temporal behavior of quantum systems is governed by unitary evolutions. 
The seminal work by Gamow on the exponential law, as well as its deriva- 
tion by Weisskopf and Wigner || are based on the assumption that a pole 
near the real axis of the complex energy plane dominates the temporal evo- 
lution of the quantum system. This assumption leads to a spectrum of the 
Breit- Wigner type || and to the Fermi Golden Rule 0. However, it is well 
known that a purely exponential decay law can neither be expected for very 
short nor for very long || times. The domain of validity of the exponential 
law is limited: the long-time power tails and the short-time quadratic behav- 
ior are unavoidable consequences of very general mathematical properties of 
the Schrodinger equation M. 

The short-time behavior ||, in particular, turns out to be very interesting, 
due to its apparently paradoxical consequences leading the so-called quan- 
tum Zeno effect. Recent theoretical and experimental work || has focussed 
on the temporal behavior of a two-level system whose Rabi oscillations, in- 
duced by an rf field, are hindered by another, "measuring" field of different 
frequency. It should be noticed, however, that the idea of making use of an 
oscillating system to test the quantum Zeno effect is at variance with the 
original proposals, based on truly unstable systems ||. For this reason, al- 
ternative schemes were recently proposed 0, PI , that do not require any 



reinterpretation of the experimental data [12 . 

The purpose of this Letter is to investigate the characteristic features of 
the short-time nonexponential region of a truly unstable system. Our atten- 
tion will be focussed on a transition of the hydrogen atom: We shall endeavor 
to give an accurate estimate of the "Zeno" time for this system. Our general 
conclusions, however, will be valid for any two-level system interacting with 
a quantum field (as far as the theory is renormalizable) . 

Let us start by outlining the main features of the problem. Let |-0o) be 
the wave function of a given quantum system at time t = 0. The evolution 
is governed by the unitary operator U(t) = exp(—iHt/h), where H is the 
Hamiltonian. The "survival" or nondecay probability at time t is the square 
modulus of the survival amplitude 

P(t) = \(4, \e-* m / h \ij )\ 2 =l-t 2 /Tt + ... (1) 
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rz 1 = ^- = \ (>o|#Vo> - (^ol^l^o) 2 ) 172 • (2) 

The short-time expansion is quadratic in t and therefore yields a vanishing 
decay rate for t — > 0. This quadratic behavior is in manifest contradiction 
with the exponential law that predicts an initial nonvanishing decay rate (the 
inverse of the lifetime). The quantity rz will be referred to as "Zeno time," 
in the present paper. 

Unfortunately, when one considers quantum field theory, things do not 
work out that easily. In the above (naive) derivation, one assumes that all 
moments of H in the state j^o) are finite and (implicitly) that |?/>o) is nor- 
malizable and belongs to the domain of definition of H If the volume of 
the box containing the system is not finite, the spectrum of the Hamiltonian 
is continuous and the Zeno time turns out to be inversely proportional to 
some power of a frequency cutoff A: rz oc 1/A a . This is a very general prop- 
erty, essentially due to the singular nature of the product of local observables 
when computed at short distances [p"4[ . 



However, if the theory is renormalizable, this divergence can be tamed by 
introducing a natural cutoff for the system. In the present paper we shall just 
concentrate our attention on such a situation: We will show that it is indeed 
possible to compute the value of rz for the 2P-1S transition of the hydrogen 
atom. The result is finite. This confirms that a quantum Zeno region is not 
a phenomenon peculiar to the quantum mechanics of finite systems; rather, 
it is present even in the more general framework of quantum field theory, at 
least for a renormalizable theory. 

We start from the total Hamiltonian (h = c = 1) 

H = H atom + Hem + -f^nt 

= 2lEi\i)(i\+22 duualpdup (3) 
i=i p Jo 

roc r , 

+ E / du J [ip p (u J )al p \l)(2\ + ip^(uj)a U)P \2)(l\\, 
p jo 

where the first term is the free Hamiltonian of a two-level atom, the second 
term the Hamiltonian of the free EM field and the third term the interac- 
tion Hamiltonian. We considered only the linear part of the interaction (in 
the so-called rotating wave approximation) and expanded it in the energy- 
angular momentum basis for photons fL5f , with J2p — S^=i Ei=-j Sa=0' 
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where |z) (i = 1, 2) are the atomic states (of energy Ei), A defines the photon 
parity P = (— 1) J+1+A , j is the total angular momentum (orbital+spin) of 
the photon, m its magnetic quantum number and 



[O'ujjmX, a lj'j'm'X'} ~ ~~ w')8jjt8 mm r8\\t. 



(4) 



The quantities <pp{uj) are the matrix elements of the interaction Hamiltonian 
between the states 



|1;1 W/J ) = \l)®\uj,j,m,X), |2;0) = |2) |0), 



(5) 



where the first ket refers to the atom and the second to the photon. We 
concentrate now our attention on the 2P-1S transition of hydrogen: |1) 
n | = = 0, mi = 0), 1 2) = \n 2 = 2, Z 2 = 1, Tn>2}- Conservation of 
angular momentum and parity ensures the validity of the selection rules 
j = 1, m = m.2, A = 1. This reduces the sum over (3 in the interaction 
Hamiltonian to the single term (3 = (1, 777,2, 1). In this case, the matrix 
elements were exactly evaluated by Moses |16|] and Seke [17 

<pp(u>) = (1,1^1^12,0) = (pp{uj)8pp 



i(xA) 



l) 



2 ^jl^mm2 8 X1 1 



(6) 



with 



A = -am e ~ 8.498 • 10 18 rad/s, 



X 



2 /2 



7T V3 



a 3 ~ 



6.435 • 10" 



(7) 



where a is the fine structure constant and m e the electron mass. A is the 
natural cutoff defining the atomic form factor and taking into account all 
retardation effects: in natural units, A = 3/2ao, where ao is the Bohr radius, 
so that wavelengths shorter than ao do not contribute significantly to the 
interaction. The physical origin of A is ascribable to the exponential behavior 
of the atomic orbitals, which fall off like exp(— r/nao) (where r is the radial 
coordinate and 77 the principal quantum number): for the 2P-1S transition, 
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the orbitals overlap like exp(— r/ao) • exp(— r/2ao) = exp(— rA). Notice that 
A is not put "by hand," like in analysis involving the dipole approximation, 
but naturally emerges from calculation [T7j . 

We assume that the system is initially (at time t = 0) in the eigenstate 
1 2, 0) of the unperturbed Hamiltonian H = -ff at0 m + ^em, whose eigenvalue 
is uj q = E 2 — Ei = \a 2 m e ~ 1.550 • 10 16 rad/s. We shall set E x = 0. It is now 
straightforward to compute the Zeno time, according to the definition (pi): 



- = (2,0|#? nt |2,0> = £ / cM(2,0|tf int |l,W>| 2 
r z - p Jo 

Inserting the values (0) of A and \ we obtain 
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rz = J— = (I) ^ -i— ~ 3 - 593 ■ 10_15s - ( 9 ) 



This is our first result. It is an estimate of the duration of the Zeno region 
for a truly unstable system. 

Observe that for hydrogen-like atoms of nuclear charge Z, the Zeno time 
scales (unfavorably) like Z~ 2 . This is because A oc Z/a and \ °c Z 2 a 3 . 

It is also worth stressing that the value of rz, due to its very structure, 
would not be modified by the presence of counter-rotating terms in the Hamil- 
tonian (|3|), whose contribution to (^) vanishes. Even the introduction of ad- 
ditional atomic levels would not modify this result, within the framework of 
the rotating wave approximation [whose validity is discussed after Eq. (27)]. 
On the other hand, a straightforward but rather lengthy calculation shows 
that the introduction of the other atomic levels, whose principal quantum 
number is z/, and of counter-rotating terms in the interaction Hamiltonian 
H[ nt yields the following expression for the Zeno time: 



r°° 1 4210 

(2,o|ff&|2,o) = x;/ cM(2,o|i4 t kW>r = — — (10) 



T /2 
' 7 



where the matrix elements are computed as in [17]. Equation ( |TDD yields a 
20% correction to the value of the Zeno time. 
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It is now interesting to look at the temporal behavior of our system at 
longer times. There is previous related work [18|, 19, 20, 21, 22] on this 
subject. The survival amplitude and its Laplace transform read 



y(t) = (2,0|e- im |2,0), 

f°° 1 

y(s) = / dte- st y(t) = {2M— ~FtM- 
Jo s + in 



(11) 



We make use of the identity 
1 



— = 1 j 1 Hm 1 

iH s + iHq s + iH in s + iH 

1 rr 1 rr 1 
1 1 

s + iH s + iH s + iH 



(12) 



and by introducing a complete orthonormal set of eigenstates of the un- 
perturbed Hamiltonian H [note that the interaction Hamiltonian H int has 
nonvanishing matrix elements only between the states (||])] we easily obtain 



Vis) 



-Q(s)y(s) 

s + iujQ 



s + iujQ + Q{s) 



Q(s) 



By inverting the transform we get 



dk\ipp(k)\' 



s + ik 



y{t) 



i 

27T2 Jb 



ds 



l 



Q(s) = —Q(sA) = -i 



dx- 



x 



[l + x 



2\4 



X — IS 



(13) 
(14) 

(15) 
(16) 



where B is the so-called Bromwich path, i.e. a vertical line at the right of all 
the singularities of y(s), and we used (§). Notice that Q and Q are self-energy 
contributions. It is straightforward to integrate ( jl6[) to get 



Q{s) 



-lbni - (88 - A8ni)s - 457ris 2 + 144^ 



+ 



96(s 2 - l) 4 
15ms 4 - 72s 5 - 3ms 6 + 16s 7 
96(s 2 - l) 4 



96 s log s 



(17) 
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The quantity Q(s) has a logaritmic branch cut extending from to —ioo, 
and no singularities on the first Riemann sheet (physical sheet). Indeed, 
the fourth order zeros of the denominator s = ±1 are also zeroes of the 
numerator and Q(±l) = ±32 7l' 7Tt ■ On the second Riemann sheet the function 
Q(s) becomes 



256 



Qui 



Q(se 



-2ni\ 



Q(s) + 2ni- 



vi : 



(18) 



where the additional term represents the discontinuity across the cut. It is 
easy to show that y(s) has a pole on the second Riemann sheet. From the 
denominator in fll^D, by expanding Qn(s) around — — 0' 



-if-0" 



we 



get a power series whose convergence radius is ?, because of the branching 
point at the origin. Therefore 



'pole 



■ij ~ xQn 



.Oi 



—I- 



xQ 

a AE 
+ 1— r— 
4 A 



- 0^ 

4 

a , 
-%— + 
4 

_ SL 

2A' 



o(x 2 ) 



(19) 
(20) 



because Qu(s) is the analytical continuation of Q(s) below the branch cut. 
By (UJ) we get 
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27t|^ )| 2 + O(x 2 ) 

+ 0( X 2 )^ 6.268 -10V 1 



i + (f) : 



P / du\ n {u)\ 2 

JO UJ — LUq 



+ O(x 2 )~0.5 X A, 



(21) 
(22) 



which are the Fermi "Golden Rule" (yielding the lifetime = 7 _1 — 1.595 • 
10 _9 s) and the second order correction to the level energy E?,. Notice that 
AE is not the Lamb shift, but only the shift of the 2P level due to its 
interaction with the ground state |23], |i"7| . Observe that for hydrogen-like 
atoms of nuclear charge Z, te oc (x^o) -1 scales like Z~ 4 , so that the ratio 
tz/te has the favorable scaling Z 2 . This might be important for experimental 
observation of the Zeno region. 
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The exponential law is readily obtained by deforming the original 
Bromwich path into a new contour C = C\ + C2, composed of a small circle 
Ci turning anticlockwise around the simple pole s po i e on the second Riemann 
sheet and a path C2 starting from —00 on the second sheet, turning around 
the branch point s = and extending back to —00 on the first sheet. We get 



y(t) = y po i e (t) +y ca t{t), 



(23) 



where 



y P oie(*) 



Ze 



'rt e -i(u> -AE)t+i£ 



Ze K = 



1 + XQ'n(Spo\e 



i + o(x) 



(24) 
(25) 



and the prime denotes derivative. Notice that x = 0(a 3 ) and ( ~ —2.02 • 
10~ 8 . As is well known, the exponential law is obtained by neglecting the 
contribution arising from the branch cut. Let us estimate the latter. From 



2/cut(*) 



— / dse sAt 
2ni Jc 2 



1 



s + if + xQ{s) 



£x -if- xQ{-&o) -if- xQn(-^o) 



where xq = 1/tA. At times t ^> A 1 (so that xq — > 0) 



ycut(*) ~ x%l 



c= 1 



xQ{o) 

5 X 



-x- 



c 



a 



l + 0( X ). 



(26) 



(27) 



(28) 



The expressions flT7j)-(|28|) are robust against the introduction of counter- 
rotating terms in the Hamiltonian (|3]), which can be shown to contribute 
only first-order corrections in x — 0(a 3 ) in ([T7]) and therefore second-order 



corrections in (EG 



and <M). 
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Summarizing, the general expressions (valid Vt > 0) for the survival am- 
plitude y(t) and survival probability P(t) = \y(t)\ 2 are respectively 

y (t) = Ze-^e~^°- AE ^ - x -?—h(t)e ir >® (29) 
CZ 

-2 X j—^e-^h(t) cos [(u - AE)t + 77 (t) - C] , (30) 
where h(t) and r)(t) are real functions satisfying 

/ .2 sin T \ 

77(0) = arctan -2— lim n(t) = 0. (31) 

\ZcosC — 1/ 

Notice the presence of an oscillatory term, in Eq. (|30|). Physically, it repre- 
sents an interesting (fully quantum mechanical) interference effect between 
the cut and the pole contribution to the survival amplitude (|29|). 
For short and long times, Eq. (|30|) yields 

P(t) -1-4 (*<t z ) (32) 
P(t) ~ ZV^ + x 2 ^-2 X ^e-i t cos[K-AE)t-C]. (33) 



(t > A 



Numerical investigation of ( p0|) shows that the "long-time" expansion is 
already valid for rather short times t > 2 • 10~ 17 s. For even shorter 
times, the system undergoes a rapid initial oscillation, of duration about 
200 ■ A -1 ~ 2.3 • 10~ 17 s, and then quickly relaxes towards the asymptotic 
expression (|3"3"D. The initial convexity of the curve is given by (|3"2"D, which 
agrees extremely well with the numerical investigation. 

The above analysis clarifies an important point: in contrast with a 
widespread, naive expectation, the short time behavior, yielding a vanishing 
decay rate, is nothing but the first of a series of oscillations, whose ampli- 
tude vanishes exponentially with time, eventually leading to a power law. 
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The asymptotic frequency of the oscillations is essentially ujq [see (pB])]: any 
correction (like our AE, or the total Lamb shift, or fine structure effects, not 
considered in our analysis) is at most of order 10~ 6 c<Jo- The transition to a 
power law occurs when the two summands in (|29|) are comparable, so that 
(uotfe-^ » x , namely for t ~ 98 lifetimes H, |0|, H & 

The above conclusions, derived for the hydrogen atom in interaction with 
the EM field, are generally valid for a renormalizable (or superrenormalizable) 
theory: Any interaction Hamiltonian of the type (|3]), which does not contain 
derivative couplings in the fields, yields similar results. One should notice, 
however, that the evaluation of the duration of the Zeno region depends on 
the frequency cutoff: in general, one expects a dependence on some inverse 
power of A |TJ]]; for example, in our case r z = 0(A _1 ), as in (g). An accurate 



estimate of A can pose in general a difficult problem. 

An interesting problem is to understand whether the initial quadratic 
behavior (|32|) is experimentally observable. This is an experimentally chal- 
lenging task, that raises interesting theoretical and experimental questions 
about the problem of state preparation. The time scales involved are very 
small, so that a sharp initial state preparation, even by modern pulsed-laser 
techniques, appears difficult. On the other hand, state preparation by means 
of indirect excitation processes, e.g. by electron or ion collision, seems more 
realistic. 

It is also worth stressing that the problem of sharply defining the initial 
moment of excitation might be circumvented: close scrutiny of Eqs. (p9|) -(p3|) 
suggests that experimental observation of the probability oscillations would 
not only provide a direct evidence of the cut contribution to the survival 
amplitude, but also an indirect, yet convincing, proof of the presence of the 
Zeno region, in the light of the discussion following (|33|). 
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